Fibrinolytic and coagulative activities of Yersinia pestis.

The outer membrane protease Pla belongs to the omptin protease family spread by horizontal gene transfer into Gram-negative bacteria that infect animals or plants. Pla has adapted to support the life style of the plague bacterium Yersinia pestis. Pla has a ß-barrel fold with 10 membrane-spanning ß strands and five surface loops, and the barrel surface contains bound lipopolysaccharide (LPS) that is critical for the conformation and the activity of Pla. The biological activity of Pla is influenced by the structure of the surface loops around the active site groove and by temperature-induced LPS modifications. Several of the putative virulence-related functions documented for Pla in vitro address control of the human hemostatic system, i.e., coagulation and fibrinolysis. Pla activates human plasminogen to the serine protease plasmin and activates the physiological plasminogen activator urokinase. Pla also inactivates the protease inhibitors alpha-2-antiplasmin and plasminogen activator inhibitor 1 (PAI-1) and prevents the activation of thrombin-activatable fibrinolysis inhibitor (TAFI). These functions enhance uncontrolled fibrinolysis which is thought to improve Y. pestis dissemination and survival in the mammalian host, and lowered fibrin(ogen) deposition has indeed been observed in mice infected with Pla-positive Y. pestis. However, Pla also inactivates an anticoagulant, the tissue factor (TF) pathway inhibitor, which should increase fibrin formation and clotting. Thus, Pla and Y. pestis have complex interactions with the hemostatic system. Y. pestis modifies its LPS upon transfer to the mammalian host and we hypothesize that the contrasting biological activities of Pla in coagulation and fibrinolysis are influenced by LPS changes during infection.

Human single-chain urokinase is activated by the omptins PgtE of Salmonella enterica and Pla of Yersinia pestis despite mutations of active site residues.

Fibrinolysis is important in cell migration and tightly regulated by specific inhibitors and activators; of the latter, urokinase (uPA) associates with enhancement of cell migration. Active uPA is formed through cleavage of the single-chain uPA (scuPA). The Salmonella enterica strain 14028R cleaved human scuPA at the peptide bond Lys158-Ile159, the site cleaved also by the physiological activator human plasmin. The cleavage led to activation of scuPA, while no cleavage or activation were detected with the mutant strain 14028R lacking the omptin protease PgtE. Complementation and expression studies confirmed the role of PgtE in scuPA activation. Similar cleavage and activation of scuPA were detected with recombinant Escherichia coli expressing the omptin genes pla from Yersinia pestis, ompT and ompP from E. coli, sopA from Shigella flexneri, and leo from Legionella pneumophila. For these omptins the activation of scuPA is the only shared function so far detected. Only poor cleavage and activation of scuPA were seen with YcoA of Y. pestis and YcoB of Yersinia pseudotuberculosis that are considered to be proteolytically inactive omptin variants. Point mutations of active site residues in Pla and PgtE had different effects on the proteolysis of plasminogen and of scuPA, indicating versatility in omptin proteolysis.

Regiospecificity and stereospecificity of human UDP-glucuronosyltransferases in the glucuronidation of estriol, 16-epiestriol, 17-epiestriol, and 13-epiestradiol.

The glucuronidation of estriol, 16-epiestriol, and 17-epiestriol by the human UDP-glucuronosyltransferases (UGTs) of subfamilies 1A, 2A, and 2B was examined. UGT1A10 is highly active in the conjugation of the 3-OH in all these estriols, whereas UGT2B7 is the most active UGT toward one of the ring D hydroxyls, the 16-OH in estriol and 16-epiestriol, but the 17-OH in 17-epiestriol. Kinetic analyses indicated that the 17-OH configuration plays a major role in the affinity of UGT2B7 for estrogens. The glucuronidation of the different estriols by the human liver and intestine microsomes reflects the activity of UGT1A10 and UGT2B7 in combination with the tissues' difference in UGT1A10 expression. The UGT1A10 mutant 1A10-F93G exhibited much higher V(max) values than UGT1A10 in estriol and 17-epiestriol glucuronidation, but a significantly lower value in 16-epiestriol glucuronidation. To this study on estriol glucuronidation we have added experiments with 13-epiestradiol, a synthetic estradiol in which the spatial arrangement of the methyl on C18 and the hydroxyl on C17 is significantly different than in other estrogens. In comparison with estradiol glucuronidation, the C13 configuration change decreases the turnover of UGTs that conjugate the 3-OH, but increases it in UGTs that primarily conjugate the 17-OH. Unexpectedly, UGT2B17 exhibited similar conjugation rates of both the 17-OH and 3-OH of 13-espiestradiol. The combined results reveal the strong preference of UGT1A10 for the 3-OH of physiologic estrogens and the equivalently strong preference of UGT2B7 and UGT2B17 for the hydroxyls on ring D of such steroid hormones.

Glutamine synthetase and glucose-6-phosphate isomerase are adhesive moonlighting proteins of Lactobacillus crispatus released by epithelial cathelicidin LL-37.

Glutamine synthetase (GS) and glucose-6-phosphate isomerase (GPI) were identified as novel adhesive moonlighting proteins of Lactobacillus crispatus ST1. Both proteins were bound onto the bacterial surface at acidic pHs, whereas a suspension of the cells to pH 8 caused their release into the buffer, a pattern previously observed with surface-bound enolase and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) of L. crispatus. The pH shift was associated with a rapid and transient increase in cell wall permeability, as measured by cell staining with propidium iodide. A gradual increase in the release of the four moonlighting proteins was also observed after the treatment of L. crispatus ST1 cells with increasing concentrations of the antimicrobial cationic peptide LL-37, which kills bacteria by disturbing membrane integrity and was here observed to increase the cell wall permeability of L. crispatus ST1. At pH 4, the fusion proteins His(6)-GS, His(6)-GPI, His(6)-enolase, and His(6)-GAPDH showed localized binding to cell division septa and poles of L. crispatus ST1 cells, whereas no binding to Lactobacillus rhamnosus GG was detected. Strain ST1 showed a pH-dependent adherence to the basement membrane preparation Matrigel. Purified His(6)-GS and His(6)-GPI proteins bound to type I collagen, and His(6)-GS also bound to laminin, and their level of binding was higher at pH 5.5 than at pH 6.5. His(6)-GS also expressed a plasminogen receptor function. The results show the strain-dependent surface association of moonlighting proteins in lactobacilli and that these proteins are released from the L. crispatus surface after cell trauma, under conditions of alkaline stress, or in the presence of the antimicrobial peptide LL-37 produced by human cells.

Molecular adaptation of a plant-bacterium outer membrane protease towards plague virulence factor Pla.

BACKGROUND: Omptins are a family of outer membrane proteases that have spread by horizontal gene transfer in Gram-negative bacteria that infect vertebrates or plants. Despite structural similarity, the molecular functions of omptins differ in a manner that reflects the life style of their host bacteria. To simulate the molecular adaptation of omptins, we applied site-specific mutagenesis to make Epo of the plant pathogenic Erwinia pyrifoliae exhibit virulence-associated functions of its close homolog, the plasminogen activator Pla of Yersinia pestis. We addressed three virulence-associated functions exhibited by Pla, i.e., proteolytic activation of plasminogen, proteolytic degradation of serine protease inhibitors, and invasion into human cells. RESULTS: Pla and Epo expressed in Escherichia coli are both functional endopeptidases and cleave human serine protease inhibitors, but Epo failed to activate plasminogen and to mediate invasion into a human endothelial-like cell line. Swapping of ten amino acid residues at two surface loops of Pla and Epo introduced plasminogen activation capacity in Epo and inactivated the function in Pla. We also compared the structure of Pla and the modeled structure of Epo to analyze the structural variations that could rationalize the different proteolytic activities. Epo-expressing bacteria managed to invade human cells only after all extramembranous residues that differ between Pla and Epo and the first transmembrane ß-strand had been changed. CONCLUSIONS: We describe molecular adaptation of a protease from an environmental setting towards a virulence factor detrimental for humans. Our results stress the evolvability of bacterial ß-barrel surface structures and the environment as a source of progenitor virulence molecules of human pathogens.

Two distinct regions in the model protein Peb1 are critical for its heterologous transport out of Escherichia coli.

BACKGROUND: Escherichia coli is frequently the first-choice host organism in expression of heterologous recombinant proteins in basic research as well as in production of commercial, therapeutic polypeptides. Especially the secretion of proteins into the culture medium of E. coli is advantageous compared to intracellular production due to the ease in recovery of the recombinant protein. Since E. coli naturally is a poor secretor of proteins, a few strategies for optimization of extracellular secretion have been described. We have previously reported efficient secretion of the diagnostically interesting model protein Peb1 of Campylobacter jejuni into the growth medium of Escherichia coli strain MKS12 (ΔfliCfliD). To generate a more detailed understanding of the molecular mechanisms behind this interesting heterologous secretion system with biotechnological implications, we here analyzed further the transport of Peb1 in the E. coli host. RESULTS: When mature Peb1 was expressed without its SecA-YEG-dependent signal sequence and without the putative signal peptidase II recognition sequence in E. coli MKS111ΔHBB lacking the flagellar secretion complex, the protein was found in the periplasm and growth medium which indicated a flagellum-independent translocation. We assessed the Peb1 secretion proficiency by an exhaustive search for transport-affecting regions using a transposition-based scanning mutagenesis strategy. Strikingly, insertion mutagenesis of only two segments, called TAR1 (residues 42 and 43) and TAR2 (residues 173 to 180), prevented Peb1 secretion individually. We confirmed the importance of TAR regions by subsequent site-specific mutagenesis and verified that the secretion deficiency of Peb1 mutants was not due to insolubility or aggregation of the proteins in the cytoplasm. We found by cell fractionation that the mutant proteins were present in the periplasm as well as in the cytoplasm of MKS12. Hence, mutagenesis of TAR regions did not affect export of Peb1 across the cytoplasmic membrane, whereas its export over the outer membrane was markedly impaired. CONCLUSIONS: We propose that the localization of the model protein Peb1 in the growth medium of E. coli is due to active secretion by a still unknown pathway of E. coli. The secretion apparently is a two-step process involving a periplasmic step and the TAR regions.

The omptins of Yersinia pestis and Salmonella enterica cleave the reactive center loop of plasminogen activator inhibitor 1.

Plasminogen activator inhibitor 1 (PAI-1) is a serine protease inhibitor (serpin) and a key molecule that regulates fibrinolysis by inactivating human plasminogen activators. Here we show that two important human pathogens, the plague bacterium Yersinia pestis and the enteropathogen Salmonella enterica serovar Typhimurium, inactivate PAI-1 by cleaving the R346-M347 bait peptide bond in the reactive center loop. No cleavage of PAI-1 was detected with Yersinia pseudotuberculosis, an oral/fecal pathogen from which Y. pestis has evolved, or with Escherichia coli. The cleavage and inactivation of PAI-1 were mediated by the outer membrane proteases plasminogen activator Pla of Y. pestis and PgtE protease of S. enterica, which belong to the omptin family of transmembrane endopeptidases identified in Gram-negative bacteria. Cleavage of PAI-1 was also detected with the omptins Epo of Erwinia pyrifoliae and Kop of Klebsiella pneumoniae, which both belong to the same omptin subfamily as Pla and PgtE, whereas no cleavage of PAI-1 was detected with omptins of Shigella flexneri or E. coli or the Yersinia chromosomal omptins, which belong to other omptin subfamilies. The results reveal a novel serpinolytic mechanism by which enterobacterial species expressing omptins of the Pla subfamily bypass normal control of host proteolysis.

A molecular model of the human UDP-glucuronosyltransferase 1A1, its membrane orientation, and the interactions between different parts of the enzyme.

The vertebrate UDP-glucuronosyltransferases (UGTs) are membrane-bound enzymes of the endoplasmic reticulum that process both endogenous and exogenous substrates. The human UGTs are well known biologically, but biophysical understanding is scarce, largely because of problems in purification. The one resolved crystal structure covers the C-terminal domain of the human UGT2B7. Here, we present a homology model of the complete monomeric human UGT1A1, the enzyme that catalyzes bilirubin glucuronidation. The enzyme can be seen as composed of four different domains: two large ones, the N- and C-terminal domains, and two small ones, the "envelope" helices and the transmembrane segment that includes the cytoplasmic tail. The hydrophobic core of the N-terminal domain and the two envelope helices that connect the large domains are shown to be structurally well conserved even among distant homologs and can thus be modeled with good certainty according to plant and bacterial structures. We consider alternative solutions for the highly variable N-terminal regions that probably contribute to substrate binding. The bilirubin binding site, known pathological mutations in UGT1A1, and other specific residues have been examined in the context of the model with regard to available experimental data. A putative orientation of the protein relative to the membrane has been derived from the location of predicted N-glycosylation sites. The model presents extensive interactions between the N- and C-terminal domains, the two envelope helices, and the membrane. Together, these interactions could allow for a concerted large-scale conformational change during catalysis.

How many and which amino acids are responsible for the large activity differences between the highly homologous UDP-glucuronosyltransferases (UGT) 1A9 and UGT1A10?

The amino acid sequences of the human UDP-glucuronosyltransferases (UGTs) 1A9 and 1A10 are 93% identical, yet there are large differences in their activity and substrate selectivity. For example, the regioselectivity in propranolol glucuronidation, the regioselectivity in dobutamine glucuronidation, and the glucuronidation rate of alpha- and beta-estradiol differ greatly between UGT1A9 and UGT1A10. To identify the residue responsible for the activity differences, we divided the N-terminal half of the two UGTs into five comparable segments by inserting four unique restriction sites at identical positions in both genes and constructing chimeras in which segments of UGT1A9 were individually replaced by the corresponding segments from UGT1A10. Activity analyses of the resulting mutants, 910A [1A10((1-83))/1A9((84-285))], 910B [1A9((1-83))/1A10((84-147))/1A9((148-285))], 910C [1A9((1-147))/1A10((148-181))/1A9((182-285))], 910D [1A9((1-181))/1A10((182-235))/1A9((236-285))], and 910E [1A9((1-235))/1A10((236-285))] indicated that more than one residue is responsible for the differences between UGT1A9 and UGT1A10. We next prepared four double chimeras, in which two of the above UGT1A9 segments were replaced simultaneously by the corresponding UGT1A10 segments. However, none of the double chimeras glucuronidated either estradiol, propranolol, or dobutamine at rates that resembled those of UGT1A10. On the other hand, studying the kinetics of 1-naphthol glucuronidation yielded more focused results, indicating that residues within segment B (84-147) contribute directly to the K(m) value for this substrate. Further mutagenesis and activity assays suggested that Phe117 of UGT1A9 participates in 1-naphthol binding. In addition, it appears that residues within segment C of the N-terminal domain, mainly at positions 152 and 169, contribute to the higher glucuronidation rates of UGT1A10.

The chemistry of the CuB site in cytochrome c oxidase and the importance of its unique His-Tyr bond.

The CuB metal center is at the core of the active site of the heme-copper oxidases, comprising a copper atom ligating three histidine residues one of which is covalently bonded to a tyrosine residue. Using quantum chemical methodology, we have studied the CuB site in several redox and ligand states proposed to be intermediates of the catalytic cycle. The importance of the His-Tyr crosslink was investigated by comparing energetics, charge, and spin distributions between systems with and without the crosslink. The His-Tyr bond was shown to decrease the proton affinity and increase the electron affinity of both Tyr-244 and the copper. A previously unnoticed internal electronic equilibrium between the copper atom and the tyrosine was observed, which seems to be coupled to the unique structure of the system. In certain states the copper and Tyr-244 compete for the unpaired electron, the localization of which is determined by the oxygenous ligand of the copper. This electronic equilibrium was found to be sensitive to the presence of a positive charge 10 A away from the center, simulating the effect of Lys-319 in the K-pathway of proton transfer. The combined results provide an explanation for why the heme-copper oxidases need two pathways of proton uptake, and why the K-pathway is active only in the second half of the reaction cycle.

A complete Neandertal mitochondrial genome sequence determined by high-throughput sequencing.

A complete mitochondrial (mt) genome sequence was reconstructed from a 38,000 year-old Neandertal individual with 8341 mtDNA sequences identified among 4.8 Gb of DNA generated from approximately 0.3 g of bone. Analysis of the assembled sequence unequivocally establishes that the Neandertal mtDNA falls outside the variation of extant human mtDNAs, and allows an estimate of the divergence date between the two mtDNA lineages of 660,000 +/- 140,000 years. Of the 13 proteins encoded in the mtDNA, subunit 2 of cytochrome c oxidase of the mitochondrial electron transport chain has experienced the largest number of amino acid substitutions in human ancestors since the separation from Neandertals. There is evidence that purifying selection in the Neandertal mtDNA was reduced compared with other primate lineages, suggesting that the effective population size of Neandertals was small.

Modeling the active-site structure of the cbb3-type oxidase from Rhodobacter sphaeroides.

The active site of the heme-copper oxidases comprises a redox-active high-spin heme and a tris-histidine copper center Cu B. Two amino acids in the close vicinity of the metals, a tyrosine and a tryptophan from helix 6, have been shown to be absolutely required for the catalytic function and should be considered part of the active site. Additionally, amino acid residues from interhelical loops strongly modify the activity. In a separate subfamily of heme-copper oxidases, the cbb 3-type oxidases, the metal centers are identical, the tyrosine is found in helix 7, but nothing is known of the corresponding tryptophan or of the involvement of the loop residues. We have observed a conserved aromatic cluster in the known oxidase structures, including the essential tryptophan and loop residues, and refined our earlier model of the cbb 3-type oxidase from Rhodobacter sphaeroides to test the feasibility of a similar structure. In the refined model, the interactions around the Delta-propionate of the high-spin heme resemble closely those seen in crystal structures of other terminal oxidases. Two alternative models (G- and C-models) that differ for the positioning of conserved tryptophans in helix 6, are presented. Molecular dynamics simulations on the catalytic subunit of the cbb 3-type oxidase model result in a conformational change of the active-site tyrosine, which may be related to different ligand-binding properties of the cbb 3-type oxidases. The relationship between sequence and functional data for defining the subfamily is discussed.

Charge parameterization of the metal centers in cytochrome c oxidase.

Reliable atomic point charges are of key importance for a correct description of the electrostatic interactions when performing classical, force field based simulations. Here, we present a systematic procedure for point charge derivation, based on quantum mechanical methodology suited for the systems at hand. A notable difference to previous procedures is to include an outer region around the actual system of interest. At the cost of increasing the system sizes, here up to 265 atoms, including the surroundings achieves near-neutrality for the systems as well as structural stability, important factors for reliable charge distributions. In addition, the common problem of converting between C--H bonds and C--C bonds at the border vanishes. We apply the procedure to the four redox-active metal centers of cytochrome c oxidase: Cu(A), haem a, haem a(3), and Cu(B). Several relevant charge and ligand states are considered. Charges for two different force fields, CHARMM and AMBER, are presented.

Dynamics of the glutamic acid 242 side chain in cytochrome c oxidase.

In many cytochrome c oxidases glutamic acid 242 is required for proton transfer to the binuclear heme a(3)/Cu(B) site, and for proton pumping. When present, the side chain of Glu-242 is orientated "down" towards the proton-transferring D-pathway in all available crystal structures. A nonpolar cavity "above" Glu-242 is empty in these structures. Yet, proton transfer from Glu-242 to the binuclear site, and for proton-pumping, is well established, and the cavity has been proposed to at least transiently contain water molecules that would mediate proton transfer. Such proton transfer has been proposed to require isomerisation of the Glu-242 side chain into an "up" position pointing towards the cavity. Here, we have explored the molecular dynamics of the protonated Glu-242 side chain. We find that the "up" position is preferred energetically when the cavity contains four water molecules, but the "down" position is favoured with less water. We conclude that the cavity might be deficient in water in the crystal structures, possibly reflecting the "resting" state of the enzyme, and that the "up/down" equilibrium of Glu-242 may be coupled to the presence of active-site water molecules produced by O(2) reduction.

Polyglutamine variation in a flowering time protein correlates with island age in a Hawaiian plant radiation.

BACKGROUND: A controversial topic in evolutionary developmental biology is whether morphological diversification in natural populations can be driven by expansions and contractions of amino acid repeats in proteins. To promote adaptation, selection on protein length variation must overcome deleterious effects of multiple correlated traits (pleiotropy). Thus far, systems that demonstrate this capacity include only ancient or artificial morphological diversifications. The Hawaiian Islands, with their linear geological sequence, present a unique environment to study recent, natural radiations. We have focused our research on the Hawaiian endemic mints (Lamiaceae), a large and diverse lineage with paradoxically low genetic variation, in order to test whether a direct relationship between coding-sequence repeat diversity and morphological change can be observed in an actively evolving system. RESULTS: Here we show that in the Hawaiian mints, extensive polyglutamine (CAG codon repeat) polymorphism within a homolog of the pleiotropic flowering time protein and abscisic acid receptor FCA tracks the natural environmental cline of the island chain, consequent with island age, across a period of 5 million years. CAG expansions, perhaps following their natural tendency to elongate, are more frequent in colonists of recently-formed, nutrient-rich islands than in their forebears on older, nutrient-poor islands. Values for several quantitative morphological variables related to reproductive investment, known from Arabidopsis fca mutant studies, weakly though positively correlate with increasing glutamine tract length. Together with protein modeling of FCA, which indicates that longer polyglutamine tracts could induce suboptimally mobile functional domains, we suggest that CAG expansions may form slightly deleterious alleles (with respect to protein function) that become fixed in founder populations. CONCLUSION: In the Hawaiian mint FCA system, we infer that contraction of slightly deleterious CAG repeats occurred because of competition for resources along the natural environmental cline of the island chain. The observed geographical structure of FCA variation and its correlation with morphologies expected from Arabidopsis mutant studies may indicate that developmental pleiotropy played a role in the diversification of the mints. This discovery is important in that it concurs with other suggestions that repetitive amino acid motifs might provide a mechanism for driving morphological evolution, and that variation at such motifs might permit rapid tuning to environmental change.

The K-pathway revisited: a computational study on cytochrome c oxidase.

Cytochrome c oxidase contains two established proton-conducting structures, the D- and K-pathways. The role of the K-pathway appears to be to conduct the first two protons to be used in water formation, which are taken up on reduction of the oxidized enzyme. Previous computational work has suggested that Lys(I)-319 is neutral over a large pH range and in various redox states. We have constructed oxidase models in different redox states using quantum-chemically derived charge parameters for the redox metal centers. The protonation behaviour of titratable sites in the two-subunit enzyme was defined by continuum electrostatics. The calculations reported here show substantial protonation of Lys(I)-319 at neutral pH once the stable X-ray crystallographic water molecule found immediately next to it is treated explicitly. The immediate structure of the Lys(I)-319 environment is independent of redox state, but the pK(a) value of this residue changes with the redox state of the binuclear heme a3/Cu(B) site whenever that change is electrically uncompensated. Lys(I)-319 is also found to interact electrostatically with the conserved residue Glu(II)-62 in subunit II. These results are discussed in relation to the role of the K-pathway in oxidase function.

Sequence analysis of the cbb3 oxidases and an atomic model for the Rhodobacter sphaeroides enzyme.

The cbb3-type oxidases are members of the heme-copper oxidase superfamily, distant by sequence comparisons, but sharing common functional characteristics. To understand the minimal common properties of the superfamily, and to learn about cbb3-type oxidases specifically, we have analyzed a wide set of heme-copper oxidase sequences and built a homology model of the catalytic subunit of the cbb3 oxidase from Rhodobacter sphaeroides. We conclude that with regard to the active site surroundings, the cbb3 oxidases greatly resemble the structurally known oxidases, while major differences are found in three segments: the additional N-terminal stretch of ca. 60 amino acids, the segment following helix 3 to the end of helix 5, and the C-terminus from helix 11 onward. The conserved core contains the active site tyrosine and also an analogue of the K-channel of proton transfer, but centered on a well-conserved histidine in the lower part of helix 7. Modeling the variant parts of the enzyme suggests that two periplasmic loops (between helices 3 and 4 and between helices 11 and 12) could interact with each other as a part of the active site structure and might have an important role in proton pumping. An analogue of the D-channel is not found, but an alternative channel might form around helix 9. A preliminary packing model of the trimeric enzyme is also presented.

Gating of proton and water transfer in the respiratory enzyme cytochrome c oxidase.

The membrane-bound enzyme cytochrome c oxidase is responsible for cell respiration in aerobic organisms and conserves free energy from O2 reduction into an electrochemical proton gradient by coupling the redox reaction to proton-pumping across the membrane. O2 reduction produces water at the bimetallic heme a3/CuB active site next to a hydrophobic cavity deep within the membrane. Water molecules in this cavity have been suggested to play an important role in the proton-pumping mechanism. Here, we show by molecular dynamics simulations that the conserved arginine/heme a3 delta-propionate ion pair provides a gate, which exhibits reversible thermal opening that is governed by the redox state and the water molecules in the cavity. An important role of this gate in the proton-pumping mechanism is supported by site-directed mutagenesis experiments. Transport of the product water out of the enzyme must be rigidly controlled to prevent water-mediated proton leaks that could compromise the proton-pumping function. Exit of product water is observed through the same arginine/propionate gate, which provides an explanation for the observed extraordinary spatial specificity of water expulsion from the enzyme.

Adaptive evolution of cytochrome c oxidase: Infrastructure for a carnivorous plant radiation.

Much recent attention in the study of adaptation of organismal form has centered on developmental regulation. As such, the highly conserved respiratory machinery of eukaryotic cells might seem an unlikely target for selection supporting novel morphologies. We demonstrate that a dramatic molecular evolutionary rate increase in subunit I of cytochrome c oxidase (COX) from an active-trapping lineage of carnivorous plants is caused by positive Darwinian selection. Bladderworts (Utricularia) trap plankton when water-immersed, negatively pressured suction bladders are triggered. The resetting of traps involves active ion transport, requiring considerable energy expenditure. As judged from the quaternary structure of bovine COX, the most profound adaptive substitutions are two contiguous cysteines absent in approximately 99.9% of databased COX I sequences from Eukaryota, Archaea, and Bacteria. This motif lies directly at the docking point of COX I helix 3 and cytochrome c, and modeling of bovine COX I suggests the possibility of an unprecedented helix-terminating disulfide bridge that could alter COX/cytochrome c dissociation kinetics. Thus, the key adaptation in Utricularia likely lies in molecular energetic changes that buttressed the mechanisms responsible for the bladderworts' radical morphological evolution. Along with evidence for COX evolution underlying expansion of the anthropoid neocortex, our findings underscore that important morphological and physiological innovations must often be accompanied by specific adaptations in proteins with basic cellular functions.

The nucleotide-binding domain of the Zn2+-transporting P-type ATPase from Escherichia coli carries a glycine motif that may be involved in binding of ATP.

In P-type ATPases, the nucleotide-binding (N) domain is located in the middle of the sequence which folds into the phosphorylation (P) domain. The N domain of ZntA, a Zn2+-translocating P-type ATPase from Escherichia coli, is approx. 13% identical with the N domain of sarcoplasmic reticulum Ca2+-ATPase. None of the Ca2+-ATPase residues involved in binding of ATP are found in ZntA. However, the sequence G503SGIEAQV in the N domain of ZntA resembles the motif GxGxxG, which forms part of the ATP-binding site in protein kinases. This motif is also found in Wilson disease protein where several disease mutations cluster in it. In the present work, we have made a set of disease mutation analogues, including the mutants G503S (Gly503-->Ser), G505R and A508F of ZntA. At low [ATP], these mutant ATPases are poorly phosphorylated. The phosphorylation defect of the mutants G503S and G505R can, however, be partially (G503S) or fully (G505R) compensated for by using a higher [ATP], suggesting that these mutations lower the affinity for ATP. In all three mutant ATPases, phosphorylation by P(i) has become less sensitive to the presence of ATP, also consistent with the proposal that the Gly503 motif plays a role in ATP binding. In order to test this hypothesis, we have modelled the N domain of ZntA using the sarcoplasmic reticulum Ca2+-ATPase structure as a template. In the model, the Gly503 motif, as well as the residues Glu470 and His475, are located in the proximity of the ATP-binding site. In conclusion, the mutagenesis data and the molecular model are consistent with the idea that the two loops carrying the residues Glu470, His475, Gly503 and Gly505 play a role in ATP binding and activation.